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Abstract 



We use path integrals to calculate hedge parameters and efficacy of hedging in a quantum field 
theory generalization of the Heath, Jarrow and Morton (HJM) term structure model which parsi- 
moniously describes the evolution of imperfectly correlated forward rates. We also calculate, within 
the model specification, the effectiveness of hedging over finite periods of time. We use empirical 
estimates for the parameters of the model to show that a low dimensional hedge portfolio is quite 
effective. 



1 Introduction 

The first interest models were spot rate models which had only one factor which implied that the prices of 
all bonds 1 were perfectly correlated. This was observed not to be the case in practice. This observation 
led Heath, Jarrow and Morton [?] to develop their famous model (heceforth called the HJM model) 
where the forward rate curve was influenced by more than one factor. This enabled bond prices to have 
imperfect correlation. However, for an iV factor HJM model, this still meant that the movements in the 
price of N bonds would determine the movements in the prices of all other bonds. This would enable 
one to hedge any instrument with N bonds within the framework of this model. However, this does not 
again seem to be the case in practice. In fact, if taken to be exact, a two factor HJM model would seem 
to imply that one can hedge a thirty year treasury bond with three month and six month bills which 
is not reasonable. Hence, there has been much interest in developing models which do not have this 
problem. One possibility is to use an infinite factor HJM model as pointed out in Cohen and Jarrow 
[?] but it is well known that estimating the parameters of even a two or three factor HJM model from 
market data is very difficult. In contrast, the estimation of parameters for different field theory models 
has been discussed in Baaquie and Srikant [?] and is seen to be more effective. 

These observations led Kennedy [?], Santa-Clara and Sornette [?] and Goldstein [?] to come up with 
random field models which allowed imperfect correlations across all the bonds. Baaquie [?] furthered 
this development by putting all these models into a field theory framework which allows for the use of a 
large body of theoretical and computational methods developed in physics to be applied to this problem. 



2 A Brief Summary of the HJM and Field Theory Models 



In the HJM model the forward rates are given by 




(1) 



where W% are independent Wiener processes. We can also write this as 




(2) 



i=l 



In this paper, we only use zero coupon bonds, hence all references to bonds are to zero coupon bonds 



where rji represent independent white noises. The action functional, is 



i=i 

We can use this action to calculate the generating functional which is 

r*2 . 



i K r 

s[w] = --J2J dt Viit) (3) 



E£i/«?*i?(*) ( 4 ) 



We now review Baaquie's field theory model presented in [?] with constant rigidity. Baaquie proposed 
that the forward rates being driven by white noise processes in (Q) be replaced by considering the forward 
rates itself to be a quantum field. To simplify notation, we write the evolution equation in terms of the 
velocity quantum field A(t,x), and which yields 

f(t,x) =f(t ,x)+ / dt , a(t',x) + Y J dt'a l (t',x)A i (t',x) (5) 

«to 2=1 

or 

df(t x) J-^ 

K dt ' = a(t, x)+J2 <7i(t, x )Mt, x ) (6) 

i=l 

The main extension to HJM is that A depends on x as well as t unlike W which only depends on t. 

While we can put in many fields A4 , it was shown in Baaquie and Srikant [?] that the extra generality 
brought into the process due to the extra argument x makes one field sufficient. Hence, in future, we 
will drop the subscript for A. 

Baaquie further proposed that the field A has the free (Gaussian) free field action functional [?] 

with Neumann boundary conditions imposed at x = t and x = t + TpR. This makes the action equivalent 
(after an integration by parts where the surface term vanishes) to 

l f°° r t+T ™ ( 1 d 2 \ 

s = ~2j t dt J t dxA & x )y--^2Q^) A & x ) ( 8 ) 

This action has the partition function 

J dt J dxdx'j{t 7 x)D(x-t,x / -t)j(t,x l )j (9) 

with 

cosh M (7> R -\9- Q'\) + cosh M (T FR -(9 + 0')) 
1 ' ' FR> ^ 2 sinh fj,T FR (10) 

= D(0', 9; T FR ) : Symmetric Function of 9, 9' 

where 9 — x — t and 9' = x' — t. We can calculate expectations and correlations using this partition 
function. Note that due to the boundary conditions imposed, the inverse of the differential operator 
D actually depends only the difference x — t. The above action represents a Gaussian random field 
with covariance structure D. In [?], a different form was found as the boundary conditions used were 
Dirichlet with the endpoints integrated over. This boundary condition is in fact equivalent to the 
Neumann condition which leads to the much simpler propagator above. In the limit Tp R — > 00 which 
we will usually take, the propagator takes the simple form ^e~^ e> cosh [x9 < where 9 > and 6*< stand for 
max(0, 9') and min(0, 9') respectively. 



When /i — > 0, this model should go over to the HJM model. This is indeed seen to be the case as 
it is seen that lim^o D(0, 6'; Tfr) = jj— • The extra factor of Tfr is irrelevant as it is due to the 
freedom we have in scaling a and D. The a we use for the different models are only comparable after 
D is normalized 2 . On normalization, the propagator for both the HJM model and field theory model in 
the limit /i — > is one showing that the two models are equivalent in this limit. 

3 Hedging 

The main aim of hedging is to reduce one's exposure to risk. There are many ways to define risk [?] . 
For bonds, the main risks are changes in interest rates and the risk of default. In this paper, we are only 
dealing with default-free bonds so that the only source of risk is the change in interest rates. 

For the purposes of this paper, we define risk to be the standard deviation or variance of final value. 
Hence, when we hedge a certain instrument, we are trying to create a portfolio of the hedged and hedging 
instruments which minimizes the overall variance of the portfolio. In the case of a N-factor HJM model, 
perfect hedging (i.e., a zero portfolio variance) is achievable once any N independent hedging instruments 
are used. However, the difficulties introduced by the infinite number of factors in the field theory models 
has resulted in their being very little literature on this important subject, a notable exception being the 
measure valued trading strategy developed in Bjork, Kabanov and Runggaldier[?]. 

In the fourth section of this paper, we will consider instantaneous hedging which is important for 
theoretical purposes. We will calculate the maximum reduction in variance for a finite number of hedging 
instruments and the hedge ratios (the amount of hedging instrument that requires to be used) that result 
in this maximization. This will show us how well the model can be approximated by a finite number of 
factors. We will then use the constant rigidity model fitted with empirical data to estimate the reduction 
in the variance of an optimally hedged portfolio as the number of hedging instruments are increased. 
We will see that a relatively small number of hedging instruments gives good results. We will also show 
that the results reduced to well known textbook ones as in Jarrow and Turnbull[?] when we go to the 
degenerate case of one-factor HJM model where all the forward rate innovations are perfectly correlated. 
We will also perform the same calculations using the propagator estimated from empirical data. 

In the third section, we will consider finite time hedging which is important in practice. This is 
because continuous hedging cannot be done in practice due to the presence of transaction costs. We will 
see how the hedging performance found in the second section changes as the time between rebalancings 
is increased. The entire analysis here is to investigate how portfolios of bonds behave in such models. 

4 Instantaneous Hedging 

In instantaneous hedging, we are considering a hedging portfolio which is rebalanced continuously in time. 
Hence, we are only considered with the instantaneous variance of the portfolio. This can be calculated 
for an arbitrary portfolio by using the fact that the covariance of the innovations in the forward rates is 
given by 

a{9)D(e,e')a(e') (11) 

in the field theory model. We will only present the hedging of zero coupon bonds in this section though 
it will be seen that the results can be easily extended to other instruments. In the first subsection, we 
will present the theoretical derivation of the hedge ratios and reduced variance for the hedging of a zero 
coupon bond with other zero coupon bonds. In the second subsection, we use the empirically fitted a 
and D(n) (from (0)) discussed in Baaquie and Srikant [?] for the constant rigidity action as well as the 
non-parametric estimate for a and D which is directly obtained from the market correlation matrix of 
the innovations in forward rates to calculate the semi-empirical reduction in variance. In the third and 
fourth subsections, we will carry out similar calculations when hedging zero coupon bonds with futures 
on zero coupon bonds. This is much more realistic in practice as hedging with futures is relatively cheap. 

2 This freedom exists since we can always make the transformation <r(£>) ~ r)(8)a(8) and D(8,8') ~ D(8, 9' ) / (^(8)^(8' )) 
without affecting any result 



4.1 Hedging bonds with other bonds 



We now consider the hedging of one bond maturing at T with N other bonds maturing at Tj, 1 < i < N . 
If one of the Ti — T, then the solution is trivial since it is the same bond. The hedge is then just 
to short the same bond giving us a zero portfolio with obviously zero variance. Since this solution is 
uninteresting, we assume that Tj ^ T \fi. The hedged portfolio 11(f) can then be represented as 

JV 
i=l 

where Aj denotes the amount of the i th bond P(t, Tj) included in the hedged portfolio. Note the value 
of bonds P(t,T) and P(t,Ti) are determined by observing their market values at time t. It is the 
instantaneous change in the portfolio value that is stochastic. Therefore, the variance of this change is 
computed to ascertain the efficacy of the hedge portfolio. 

We first consider the variance of the value of an individual bond in the field theory model. The 
definition P(t, T) = exp (— J, dxf(t, x)) for zero coupon bond prices implies that 

^|-y = f(t,t)dt-j T 'dBdf&O) 



and E 



dP{t,T) 
P(t.t) 



/•(/)-/' d0a(6)- [ d0a(6)A(t,9)) dt 



= (r(t) - J T * d6a(6)) dt since E[A(t, 9)} = 0. Therefore 



dP(t,T) 
P(t,T) 



~dP(t,T)~ 


= -dt [ 


[ PfaT) _ 


Jo 



i-T-t 

d0a{9)A(t,e) (12) 



Squaring this expression and invoking the result that E[A(t, 0)A(t, 6')] = 6(0)D(6, 0'\T FR ) = D(9 ' 8 '^ fr) 
results in the instantaneous bond price variance 



j-T-t pT-t 

Var[dP(t,T)] = dtP 2 (t,T) d9 dd'a{6)D(9,9';T FR )a(9') (13) 

Jo Jo 

As an intermediate step, the instantaneous variance of a bond portfolio is considered. For a portfolio 
of bonds, II(t) = 2i=i AjP(f, Tj), the following results follow directly 

N rTi-t 

dll(t) - E[dn(t)] = -dt V AiP(i, Ti) I d9a{9)A(t, 9) (14) 

and 

N N rTi-t pTj-t 

Var[<m(t)]=dt^y2A i A j P(t > Ti)P(t ) T j ) d9 d9'a{9)D(9,9';T FR )a(9') (15) 

i=i j=i Jo Jo 

The (residual) variance of the hedged portfolio 

JV 

n(t) = P(*,T) + 5^A i P(t,T i ) (16) 



i=i 



may now be computed in a straightforward manner. For notational simplicity, the bonds P(f, Tj) (being 
used to hedge the original bond) and P(t,T) are denoted Pf and P respectively. Equation (15) implies 



the hedged portfolio's variance equals the final result shown below 

rX-t r T-t 



pi —i pi —z 

P 2 d9 d9'a{9)a(9')D(9,9';T FR ) 

Jo Jo 

+2PJ2A 1 P 1 J^ d9 jf ' d9 l a(9)a{9')D(9 1 9' ] T FR ) 

N N „Ti-t fTj-t 

+ VyA,A/ 1 P J d9 d9'*(6)D(6,0';T FR )a(9') 

i=i j=i J ° J ° 

Note that the residual variance depends on the correlation structure of the innovation in forward rates 
described by the propagator D. Ultimately, the effectiveness of the hedged portfolio is an empirical 
question since perfect hedging is not possible without shorting the original bond. This empirical question 
is addressed in the next subsection where the propagator calibrated to market data is used to calculate 
the effectiveness. Minimizing the residual variance in equation ( |l7|) with respect to the hedge parameters 
Aj is an application of standard calculus. We introduce the following notation for simplicity. 

/■T—t i-Ti-t 

U = PP l / d9 d9'<j(6)a{9')D{e,6';T FR ) 
Jo Jo 

My = PiPj d9 / 3 d9'a{9)a(9')D{9,9';T FR ) 

Jo Jo 

Li is the covariance between the innovations in the hedged bond and the ith hedging bond and Mjj is 
the covariance between the innovations of the ith and jth hedging bond. 



The above definitions allow the residual variance in equation (|17| ) to be succinctly expressed as 

rT-t pT-t N N N 

P 2 d9 d9'a(9)<j(9 / )D(9,9';T FR ) + 2Y / ^+Y,Y, A * A 3 M v ( 18 ) 
Jo Jo i=i i=i j=i 

The hedge parameters in the field theory model can now be evaluated using basic calculus and linear 
algebra to obtain 

N 

and represent the optimal amounts of P(t, T{) to include in the hedge portfolio when hedging P(t, T). 
Putting the result into ([l7|), we see that the variance of the hedged portfolio equals 

T — t T — t N N 

V = P 2 [ d9 [ d9 l a{9)a{9')D{9,9'-T FR )-y2y2L l Mr j 1 L J (20) 
Jo Jo i=l j=1 

which declines monotonically as N increases. 

The residual variance enables the effectiveness of the hedged portfolio to be evaluated. Therefore, 
this result is the basis for studying the impact of including different bonds in the hedged portfolio as 
illustrated in the next subsection. For N = 1, the hedge parameter reduces to 

P ( tf- t d9j^- t d9'a(9)a(9')D(9,9';T FR ) \ 
p i\J Tl ~ t d9fi 1 - t d9>a(9)a(9>)D(9,9>- 1 T FR )J 

To obtain the HJM limit, we let the propagator equal one. The hedge parameter in equation (^) then 
reduces to 



P \ ^~ t d9^- t d9'a{9)a{9') \ P ( d9a{9) \ 
P ' I (j^- t d9a(9)) 2 ) ~ Pl \^- t ma{9)) 



(22) 



The popular exponential volatility function 3 a (6) = ae xe allows a comparison between our field theory 
solutions and previous research. Under the assumption of exponential volatility, equation (22) becomes 



Ai 



P_ f 1 - e -*( T -*) 
Pi I l-e-HTi-t) 



(23) 



Equation (|23j) coincides with the ratio of hedge parameters found as equation 16.13 of Jarrow and 
Turnbull [?]. In terms of their notation 



Ai 



P(t,T) f X(t,T) 
P(i,l\) \X{t,T{) 



(24) 



For emphasis, the following equation holds in a one factor HJM model 4 

0[P(t,T) + AiP(t,Ti)] 



dr(t) 



= 



(25) 



which is verified using equation ( p4| ) and results found on pages 494-495 of Jarrow and Turnbull [?] 
d[P(t,T)+A 1 P(t,T 1 ) 



dr(t) 



= -P(t,T)A:(t,T)-AiP(t,T 1 )A'(t,T 1 ) 
= -P(t,T)X(t,T)+P(t,T)X(t,T) = 



When T\ = T, the hedge parameter equals minus one. Economically, this fact states that the best 
strategy to hedge a bond is to short a bond of the same maturity. This trivial approach reduces the 
residual variance in equation (fl8|) to zero as Ai = — 1 and P = Pi implies L\ = M\\. Empirical results 
for nontrivial hedging strategies are found in the next subsection where the calibrated propagator is used. 



Sigma estimated from data 




Time to maturity/year 

Figure 1: Implied volatility function (unnormalized) for constant rigidity model using market data 



4.2 Semi-empirical results : Constant Rigidity model 

The empirical estimation of parameters for the field theory model was explained in detail in Baaquie and 
Srikant [?]. For this subsection, we use the function a estimated for the constant rigidity model from 
market data. The function a is plotted in figure [j]. 

This approach preserves the closed form solutions for hedge parameters and futures contracts illus- 
trated in the previous subsection. However, the original finite factor HJM model cannot accommodate 
an empirically determined propagator since it is automatically fixed once the HJM volatility functions 
are specified. Later in this subsection, we will see how the empirical propagator modifies the results 

3 This volatility function is commonly used as it lets the spot rate r(t) follow a Markov process. See [?]. 
4 Note that this result depends on the fact that the spot rate r(t) is Markovian and therefore only applies to either a 
constant or exponential volatility function. 



Correlation between different maturities 




Figure 2: Propagator Implied by the constant rigidity field theory model with fj, = 0.06/year 



Residual variance when hedging a 5 year zero coupon bond 
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Figure 3: Residual variance for five year bond versus bond maturity used to hedge 



of this subsection. The implied propagator for the empirically fitted value of 0.06 /yr for fi is shown in 
figure H 

The reduction in variance achievable by hedging a five year bond with other bonds is the focus of this 
subsection. We take the current forward rate curve to be flat and equal to 5% throughout. The initial 
forward rate curve does not affect any of the qualitative results. The results can also be easily extended to 
other bonds. The residual variances for one and two bond hedged portfolios are shown in figures |and§ 
The calculation of the integrals involved was done using simple trapezoidal integration as the data is not 
exceptionally accurate in the first place. Secondly and more importantly, the errors involved will largely 
cancel themselves out, hence the difference in the variances is still quite accurate. For example, in figures 
|3| and [|, we can see that in the case of perfect hedging, we get exactly zero residual variance which shows 
that the errors tend to cancel. The parabolic nature of the residual variance is because /i is constant. A 
more complicated function would produce residual variances that do not deviate monotonically as the 
maturity of the underlying and the hedge portfolio increases although the graphs appeal to our economic 
intuition which suggests that correlation between forward rates decreases monotonically as the distance 
between them increases as shown in figure ^[ Observe that the residual variance drops to zero when the 
same bond is used to hedge itself, eliminating the original position in the process. The corresponding 
hedge ratios are shown in figure |^. 

It is also interesting to note that hedging by two bonds, even very closely spaced ones, seems to be 
bring significant additional benefits. This can be seen in figure ^ where the diagonal 9 = 9' represents 
hedging by one bond. The residual variance there is higher than the nearby points in a discontinuous 
manner. 

We now present the results for the actual propagator found from the data which is graphed in figure pi 
The residual variance when a five year bond is hedged with one and two bonds bond is shown in figureslq 
and |^. We can see from figure || that, when the market propagator is used, the advantage of using more 



Residual variance when hedging a 5 year zero coupon bond 
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Figure 4: Residual variance for five year bond versus two bond maturities used to hedge 

Hedge ratio when hedging a five year zero coupon bonds 
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Figure 5: Hedge ratios for five year bond 



than one bond to hedge is significantly higher. This is because of the nature of the correlation structure 
in figure §. We see that the correlations of innovations of nearby forward rates of higher maturity is 
significantly higher in the market propagator, making hedging with more than one bond more useful. 
This is even more pronounced when hedging a short maturity bond with longer maturity ones. We can 
see this from figure [w] which shows the residual variance when a one year bond is hedged with two bonds 
where the calculation is done using the empirical propagator. The effect of this higher correlation among 
forward rates of higher maturity can also be seen in figure ^ where the residual variance rises much more 
slowly when the hedging bonds are chosen to be of higher maturity than the hedged bond. 



4.3 Hedging with futures 



We can carry out an analysis very similar to subsection 4.1 to find the optimal hedge ratios when hedging 
a bond with futures contracts on the same or other bonds. In this case, there is no trivial solution to the 
hedging problem as when bonds were hedged with other bonds. Further, since this method of hedging is 
much more practical in reality, the results will be more interesting. Proceeding as in subsection 4.1, we 
compute the appropriate hedge parameters for futures contracts. The futures price ^(tyt^jT) in terms 
of the forward price pj^'f^ = e~ - * dOffaO) an( ^ deterministic quantity Qjr(t,t*,T) which is given 

by [?] 

n^(t ,U,T) = - V / dt d6<ji(t,6) de'cr t (t,e') (26) 

■ =1 Jt Jo Jt,-t 

The dynamics of the futures price dJ-(t,t*,T) is thus given by 



dF(t,U,T) _ dn ^ t )4)T )_ f T ' d 6df{t,G) (27) 

Jt,~t 



Observed Correlation 




Figure 6: The propagator implied by the market data 



Volatility implied by the market 
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Figure 7: The volatility implied by the market data when using the empirical propagator 



d Ht^T)-E [dH tM_ 2 n = _ dt r 1 d9a[e)AM (28) 



which implies 

dF{t, i*, 

F(t,U,T) Jtm _ t 
Squaring both sides as before leads to the instantaneous variance of the futures price 

j-T-t pT-t 

Var[dF(t,U,T)]= dtF 2 (t,U,T) d6 d6'a{p)D(e,6')a{9') (29) 

Jt,-t Jt,-t 



Let Ti denote the futures price J-(t,t*,Ti) of a contract expiring at time t* on a zero coupon bond 
maturing at time Tj. The hedged portfolio in terms of the futures contract is given by 



N 



n(t) = P + (30) 



where T% represent observed market prices. For notational simplicity, define the following terms 

/•Ti-t r T-t 



Li = PTi / d0 de'a{6)D{0,e';T FR )(j{8') 
Ju-t Jo 

M i3 = T,T 3 I dO I d6'a{6)D{6,6' ;T FIi )cr{6') 
Jt.—t Ju—t 



The hedge parameters and the residual variance when futures contracts are used as the underlying 
hedging instruments have identical expressions to those in ( |l9|) and (p0|) but are based on the new 
definitions of Lj and My above. Computations parallel those in section 4.1. 



Residual variance when a 5 yr ZCP is hedged with a T yr ZCP 




Figure 8: The residual variance when a five year bond is hedged with one bond 

To explicitly state the results, the hedge parameters for a futures contract that expires at time t* on 
a zero coupon bond that matures at time Tj equals 

N 

while the variance of the hedged portfolio equals 



pT-t pT-t 

V = P 2 d9 

Jo Jo 



N N 



dO'a(6)a(9')D(e, 9'; T FR ) - X) £ LiM ij % 

i=i j=i 



for Li and M^- as defined above. 



4.4 Semi-empirical results for hedging with futures 

We first present results for the propagator fitted for the constant rigidity model as for the bonds. The 
initial forward rate curve is again taken to be flat and equal to 5%. We also fix the expiry of the futures 
contracts to be at one year from the present. This is a long enough time to clearly show the effect of the 
expiry time as well as short enough to make practical sense as long term futures contracts are illiquid 
and unsuitable for hedging purposes. 

The calculations were done using simple trapezoidal integration as explained previously. This is 
sufficient for our purposes as the fitted values for a and D shown in figures ^ and ^| are reasonably 
but not exceptionally accurate and we are more interested in the qualitative behaviour of the residual 
variance and hedge parameters. 



Residual variance when a 1 yr ZCP is optimally hedged with a T and T' yr ZCP 
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Figure 9: The residual variance when a five year bond is hedged with two bonds 



Residual variance when a 1 yr ZCP is optimally hedged with a T and T yr ZCP 




Figure 10: The residual variance when a one year bond is hedged with two bonds 



Residual variance when hedging a 5 year ZCP with one futures contract 
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Figure 11: Residual variance for a five year bond hedged with a one year futures contract on a T maturity 
bond 



The residual variance achieved when hedging a five year bond with one futures contract is shown in 
figure [H]. The optimal hedge ratios and the resulting residual variances when hedging with two and three 
futures are shown in table [l| These were obtained by systematically tabulating all possible combinations 
of bonds with intervals of three months in the maturity direction, tabulating the residual variance for 
each and finding the best combination. 

Firstly we note that the hedging is very effective even when one futures contract is used reducing the 
variance by a factor of over three hundred. Secondly, we also note that the most effective hedging is not 
obtained by shorting the futures corresponding to the same bond but one with a slightly lower maturity. 
This is due to the correlation structure of the forward rates. However, when two futures contracts are 
used, we see that one of the optimal contracts is the future on the same bond as well as a very short 
maturity futures which is probably due to the short end of the forward rate curve which does influence 
the bond but not the futures. Since the shortest maturity futures contract is probably likely to have the 
highest correlation with this part of the forward rate curve, it seems reasonable to select this futures 
contract to balance the effect on the bond from this part of the forward rate curve. This is indeed the 
case as seen in table [j]. We also note that there is very little extra improvement as we use more than 
two futures. 

We now present the same results using the empirical propagator directly. The residual variance when 



one futures contract is used for hedging is shown in figure 12. The optimal hedging futures, hedge ratios 
and residual variances are shown in table ^. We see that for the actual propagator, the optimal hedging 
futures are even farther from the actual underlying bond when compared to the optimal values using the 
fitted propagator. 



Number 


Futures Contracts (Hedge Ratio) 


Residual 
ance 


Vari- 





none 


1.82 x 10" 


-3 


1 


4.5 years (-1.288) 


5.29 x 10" 


-6 


2 


5 years (-0.9347), 1.25 years (-2.72497) 


1.58 x 10" 


-6 


3 


5 years (-0.95875), 1.5 years (1.45535), 1.25 years (-5.35547) 


1.44 x 10" 


-6 



Table 1: Residual variance and hedge ratios for a five year bond hedged with one year futures contracts. 



Residual variance when a 5 yr ZCP is optimally hedged with a T yr ZCP future 
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Figure 12: Residual variance when a five year bond is hedged with a one year futures contract on a T 
maturity bond 



5 Finite time hedging 

The case of finite time hedging is considerably more complicated. We will only do the hedging of bonds 
with other bonds as the calculations for minimizing variance can be done exactly. We will not do hedging 
of bonds with futures even though this can also be solved exactly for minimizing the variance as it does 
not add much extra insight for finite time. To see this, consider hedging with a futures contract on a 
zero coupon bond of duration T that matures at the same as the hedging horizon. This gives exactly 
the same result as hedging with a bond of the same maturity T. Therefore, we gain nothing by carrying 
out that calculation. 

The following calculation proceeds efficiently because of the use of path integral techniques which 
are very useful for such problems. To be able to optimally hedge bonds with other bonds in the sense 
of having a minimal residual variance, we need to the covariance between the final values of bonds of 
different maturities. To calculate this covariance, we will first find the joint probability density function 
for N bonds at the hedging horizon. Let us denote the initial time by 0, the hedging horizon by t\ and 
the maturities of the bonds by Tj . Making use of (0) we obtain the joint distribution of the quantities 
Gi = J t * dx(f(ti 1 x) — /(0,x)) 5 which represent the logarithms of the ratios of final value of the bonds 
to the value of their forward prices for the final time at the initial time. In other words 

= /P (fl ,r)P(*o, fl )\ = / P(t u T) \ 

\ P(t Q ,T) J \F(to,h,T)J V ' 

The calculation proceeds as follows 

N 



f[S([ 3 dx(f(t 1 ,x)-f(Q,x))-G 3 )) 

3 = 1 Jt * 

f ( * ( /•*! r T i r*i r T s 

— I dpjDAexp iy Pj \ dt / dxa(t,x) + dt dxa(t,x)A(t,x) — Gj 



(32) 



5 Due to the definition of this quantity, it is easier to carry out the calculations for the finite case with the x rather than 
the 9 variable, hence we use this variable in this section 



Number 


Futures Contracts (Hedge Ratio) 


Residual 
ance 


Vari- 





None 


1.74 x 10" 


-3 


1 


4.25 years (-0.984) 


6.34 x 10" 


-6 


2 


1.25 years (-3.84577), 5.5 years (-0.76005) 


2.26 x 10" 


-6 


3 


1.25 years (-8.60248), 1.5 years (2.84177), 5.25 (-0.85915) 


1.95 x 10" 


-6 



Table 2: Residual variance and hedge ratios for a five year bond hedged with one year futures contracts. 



which, on applying (ph becomes 



N N 



dpj exp - - p i pk / dt 



j=l k=l 



dx I dx'a(t,x)D(x — t,x' — t)a(t,x') 



N 



/ dt dxa(t, x) — Gj 
Jo Jti 



(33) 



Performing the Gaussian integrations, we obtain the joint probability distribution given by 

/ N N \ 

(27r)-"/ 2 (det J B)- 1 /2 ex p _± £ £(G j - mj )B^{G k - m k ) 

\ 3=1 fe=l / 

where B is the matrix whose elements Bij are given by 



(34) 



B^ - 



and rrii is given by 



dt j dx I dx'a{t,x)D(x — t,x' — i)cr(t,x') 

Jtx Jti 



/ dt dxa(t,x 
Jo Jti 



(35) 



(36) 



Hence, the quantities Gi follow a multivariate Gaussian distribution with covariance matrix B^ and 
mean m,. 

Having found the joint distribution of Gi, we can find the covariance of the final bond prices by 
tabulating the expectations of each of the bonds and the expectation of their products. The final bond 
price is given by P(ti,Ti) = F(0,ti,Ti)e * in terms of Gi. Hence, the expectation of this quantity is 
given by 

(2n)- N / 2 (detB)- 1 / 2 J dG,F(0, h, Ti)e Gi exp f-~(G - m) T B^(G - m)j (37) 

which gives ^(O, t\, T{) as it must since the expectation of the future bond price is the futures price. The 
expectation of the products of the prices of two bonds (P(ti,Ti)P(ti,Tj)) is given by 



(38) 



(2 7 r)- Ar / 2 (detB)- 1 / 2 J dG l dG J F(Q,t 1 ,T,)F(Q,t 1 ,T J )e G ^ +G ^ exp (—^H T B~ l H 
where H stands for the vector G — m. On evaluation, this gives the result 

JX0,*i,Ti);F(0,ii,T.,)exp ( / dt j dx [ ' dx'cr(t,x)D(x - t,x' - t)a(t,x')\ (39) 



We now consider the behaviour of the portfolio 



/v 



P(t,T) + J2^P(t,T t ) 



i=l 



(40) 



The covariance between the prices P{t\,Ti) and P(U,Tj) is given by 



exp i^j dt dx £ dx'a(t,x)D(x - t,x' - t)a(t,x')^j - ij 



(41) 



and the covariance between the hedged bond of maturity T and the hedged bonds is given by 
U =^(0,ti,T)^(0,ti,Ti) 

exp M dt £ dx dx'a(t, x)D(x —t,x' — t)a(t, x') — 1 



(42) 



and the minimization of the residual variance of the hedged portfolio proceeds exactly as in the first 
section. The hedge ratios are found to be given by 



and the minimized variance is again 



A = L T M~ X 



Nzx\P(t,T)\- L T M~ X L 



(43) 



(44) 



It is not too difficult to see that both M and L reduce to the results in the first section if — s- (with 
the covariances being scaled by t\, of course). 

One very interesting difference between the instantaneous hedging and finite time hedging is that 
the result depends on the value of a. In the calculation above, we used the risk-neutral a obtained 
for the money market numeraire. However, the market does not follow the risk-neutral measure and 
it would be better to use a value for a estimated for the market for any practical use of this method. 
This difference is expected since in the very short term only the stochastic term dominates making the 
drift inconsequential. This, of course, is not the case for the finite time case where the drift becomes 
important (it is not difficult to see that the importance of the drift grows with the time horizon) . 



Hedging a five year ZCP using one other bond over a one year interval 



Maturity of hedging bond/yr 



Figure 13: Residual variance when a five year bond is hedged with one other bond (best fit of the constant 
rigidity field theory model) with a time horizon of one year 



6 Semi-empirical Results for Finite Time Hedging 

We now present the empirical results for hedging of a bond with other bonds for both the best fit for 
the constant rigidity field theory model and the fully empirical propagator. The calculation of L and M 
were again carried out using simple trapezoidal integration and a was assumed to be purely a function 
of = x — t so that all the integrals over x were replaced by integrals over 9. The bond to be hedged 
was chosen to be the five year zero coupon bond and the time horizon t\ was chosen to be one year. 







Hedge ratio for hedging a five year ZCP using other bonds over a one year interval 
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Figure 14: Hedge ratio when a five year bond is hedged with one other bond (best fit of the constant 
rigidity field theory model) with a time horizon of one year 



Figure 15: Residual variance when a five year bond is hedged with two other bonds (best fit of the 
constant rigidity field theory model) with a time horizon of one year 

The results for the best fit of the constant rigidity field theory model (see figures p] and ||) for the 
residual variance and hedge ratio for hedging with one bond are shown in figures [K5] and j^J. The residual 
variance for hedging with two bonds is shown in figure 

The results for the fully empirical quadratic fit (see figures || and ^) for the residual variance and 
hedge ratio for hedging with one bond are shown in figures [| and 0. The residual variance for hedging 
with two bonds is shown in figure 

One interesting result is that the actual residual variance after hedging over a finite time horizon 
is lesser than naively extrapolating from the infinitesimal hedging result. This seems to be due to the 
shrinking nature of the domain as the contribution to the variance of the bonds reduces as the time 
horizon increases. This is very clear if the maturity of the bond is close to the hedging horizon as the 
volatility of bonds reduces quickly as the time to maturity approaches. Apart from this reduction, the 
results look very similar to the infinitesimal case. This is probably due to the fact that the volatility 
is quite small so the nonlinear effects in the covariance matrix ( fi"l| ) are not apparent. If very long time 
horizons (ten years or more) and long term bonds are considered, the results will probably be quite 
different. We see by comparing figures ^ and ^ that the better improvement in using more than one 
bond to hedge when the empirical rather than the field theory model propagator is used is seen to be 
true in the finite time case as well. 



Hedging a five year ZCP using two other bonds over a one year interval 
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0.00025 



Hedging a five year ZCP using one other bond over a one year interval 
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Figure 16: Residual variance when a five year bond is hedged with one other bond (best empirical fit) 
with a time horizon of one year 

Hedge ratio for optimal hedging of a five year ZCP using one other bond over a one year interval 
I 1 1 1 1 1 1 1 
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Figure 17: Hedge ratio when a five year bond is hedged with one other bond (best empirical fit) with a 
time horizon of one year 

7 Conclusion 

We have shown that the field theory model offers techniques to calculate hedge parameters for fixed 
income derivatives and provides a framework to answer questions concerning the number and maturity 
of bonds to include in a hedge portfolio. We have also seen how the field theory model can be used 
to estimate hedge parameters for finite time as well which is useful in practice. We have used the field 
theory model calibrated to market data to show that a low dimensional basis provides a reasonably 
good approximation within the framework of this model. This shows that field theory models address 
the theoretical dilemmas of finite factor term structure models and offer a practical alternative to finite 
factor models. 
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Hedging a five year ZCP using two other bonds over a one year interval 
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Figure 18: Residual variance when a five year bond is hedged with two other bonds (best empirical 
with a time horizon of one year 



